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Abstract — The discrete cosine transform (DCT), introduced by 
Ahmed, Natarajan and Rao, has been used in many applications 
of digital signal processing, data compression and information 
hiding. There are four types of the discrete cosine transform. 
In simulating the discrete cosine transform, we propose a 
generalized discrete cosine transform with three parameters, and 
prove its orthogonality for some new cases. A new type of discrete 
cosine transform is proposed and its orthogonality is proved. 
Finally, we propose a generalized discrete W transform with three 
parameters, and prove its orthogonality for some new cases. 
Keywords: Discrete Fourier transform, discrete sine transform, 
discrete cosine transform, discrete W transform 



I. Introduction 

Discrete Fourier transform has been an important tool in 
many applications of digital signal processing, image process- 
ing and information hiding. The appearance of fast fourier 
transform (FFT) has greatly promoted the rapid development 
of the subjects above. Ahmed, Natarajan, and Rao (1974) 
proposed discrete cosine transform defined on real number 
field, it can be called DCT-II-E or DCT-III-E (Wang and 
Hunt,1983). Jain (1974)proposed discrete fourier transform 
DCT-IV-E, and Wang and Hunt (1983)proposed discrete co- 
sine transform DCT-I-E. The discrete cosine transform has 
been used in frequency spectrum analysis, data compression, 
convolution computation and information hiding. Its theory 
and algorithms have received much attention for the last two 
decades (Wang and Hunt, 1983, Wang, 1984, August,2004, 
Kunz,2008). 

It is demonstrated that the performance of discrete co- 
sine transform can well approximate to ideal K-L trans- 
form (Karhunen-Loeve Transform) (Ahmed,Natarajan and 
Rao, 1974). K-L transform was proposed to dealing with a 
class of extensive stochastic image. After the image being 
transformed with K-L transform, the image restored from the 
result is the best approximation to the original image in the 
statistical sense. Moreover, for the common data model of 
Markov process, when the correlation coefficient r = 1, K- 
L transform is degraded to the classic DCT transform. In 
fact, Real-world images are neither stationary nor Markovian. 
They have different textures and structures, important image 
structures like edges, arris and lines extend over large dis- 
tances in the image (Kunz,2008). Therefore, different types of 
transformation are desirable to meet the different applications. 

Actually, the discrete cosine transform (DCT-III-E) can be 
generalized to the unified form with parameters p, q and r, as 



follows: 



JV-l 



X(k) = a(k) x(n) cos 



n=0 



k(Aqn + r)pn 
2N ' 
fc = 0,l,--- ,N-1. 



There are many new meaningful transforms, such as when 
(p,q,r) — (1, 1, 1), the new transform is orthogonal. Gener- 
ally, when gcd(p<7, N) — 1, gcd(pr, 2) = 1, where p, q and r 
are positive integers, the new transform is orthogonal. 

In a similar way, we generalize discrete cosine transforms 
DCT-II-E and DCT-IV-E. Furthermore, a new type of discrete 
cosine transform, a new type of discrete sine transform and a 
new type of discrete sine-cosine transform are proposed, and 
their orthogonality are proved. 

The discrete W transform (DWT), introduced by Wang Z, 
has been used in many applications of digital signal process- 
ing, data compression and information hiding. The unified 
form of discrete W transform has two parameters. There are 
four useful types of the discrete W transform. In simulating 
the discrete W transform, we propose a generalized discrete W 
transform with three parameters, and prove its orthogonality 
for some new cases. 

II. Generalized discrete cosine transform 

Let {x(n); n = 0, 1, 2, • • • , N} be a vector of real numbers. 
The definitions of four common types of discrete cosine 
transform (Wang and Hunt,1983) are given as follows: 

DCT-I-E: 



X(k) 
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DCT-II-E: 



JV-l 



X(k) — a(n)x(n) cos 



(2k + l)mr 



2N 

fe = 0,l,-.< ,N- 



where a(n) = 



DCT-III-E: 



n = 



N 

w else 



JV— 1 



X(k)=a(k) J2 x ( n ) 

cos 



71=0 



k(2n+ 1)tt 
2TV ' 



fc = 0,l,--- ,7V-1, 



where a(fc) = 
DCT-IV-E: 



f else 



n cos ■ 



(2fe + l)(2n + l)7r 
2TV ' 
fc = 0,l,--- ,7V-1. 



Actually, the discrete cosine transform (DCT-III-E) can be 
generalized to the following unified form with parameters p, 
q and r: 

X(k) = a(k) N f X (n)cos k ^ n 2 + r ^ , (1) 

71=0 

fc = 0,l,--- ,N-l. 

We now prove that when gcd(pq, TV) = 1, gcd (pr, 2) = 
1, where p, q and r are positive integers, transform (0 is 
orthogonal. 

Transform © can be written into matrix form: 

X(N) = C(TV) ■ x(N) 
where X(N), x(N) are column vectors of length TV, 



C < N > = Jw 



V2 
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2r P 7T 2(4g + r) P 7i 

2N 2N 



»7(4g + r-) 
2iV 
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2(8q + r)p J 
2JV 



i?(8g + r) 
2JV 



V2 

B lM«_^±Ll£i 

2[4 9 (JV-l) + r] pJ 



3 l,j(iV-l) + r] 
2N 



where rj = (TV — l)pn 

To prove that the new transform (1) is orthogonal, or 
equivalently, to prove that C(N) is an orthogonal matrix. For 

the sake of simplicity, the coefficients ^J~^ and ^= is omitted. 

Let < fci < fc 2 < TV, the inner product of the kith row 
and fc 2 th row of C(N) is that 

kirpn k 2 rpir 
cos — : cos 



27V 27V 
ki (Aq + r)pn fc 2 (4g + r)pir 



+ 



2TV 2TV 
fci [4<?(TV - 1) + r]p?r fc 2 [4g(TV - 1) + r]pn 

~\ COS COS ~~ , , 

2iV 2TV 
1 (fci - k 2 )rp-n (ki + k 2 )rpn 

771 COS — : h COS — 

2 1 27V 27V 

(ki - k 2 )(4q + r)pn (ki + k 2 )(4q + r)pn 

+ cos h cos 

27V 27V 

(h - fc 2 )[4q(TV- 1) +r]pn 
+ --- + cos — 

(fci + fc 2 )[4q(TV- 1) +r] P 7T, 
+ C ° S 2TV } 



where the above equality follows by 

cos A cos B = - [cos(A - B) + cos(A + B)] 

Since e lx = cos x + i sin x, we consider the real part of the 
following complex number: 



- krp-K ■ k(4q + r) P TV 

e 2 « + e 2 « H 1 

z l ^m- i 2k 3£2L N_ 

~ ,; 2k qP 7r U ~ e ™ J 

1 — e « 



, fc[4 g (JV-l) + r]pn 

e 2JV 



Note that gcd(pg, TV) = 1 and gcd(pr, 2) = 1. When < 

- 2fcpq7T 

k < 27V, if 1 - e l ~N— = 0, then k = TV, the real part of 



■ fcrpTr - fc(4g+T-)p7T 

2JV -)- g £ 2N 



is zero. 

If 1 - e 1 ^ 1 ^ 0, then 

■ krp-K ■ fc(4g + r-)pir 

e 2JV -|_ g J 2JV 



.. fc[4g(jV-l) + r]pn 

• e 2jv 



■ t[4i(»-l) + r]p, 

e 2JV =0 



From equality (0, when < fci < k 2 < TV, the inner 
product of the fcith row and fc 2 th row of C(TV) is zero. We 
know that C(7V) is an orthogonal matrix. 

Let < fci < TV, the inner product of the fcith row and 
itself of C(TV) is that 

2 kirpm 2 fci(4g + rW 2 
— ( (cos ) + cos + • • ■ 

7V U 27V ' 1 27V J 
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Note that 



• krpn - fc(4<; + r)pir 

e « + e " H 1 



. fc[4 g (JV-l) + r]pn 



1 — e « 

and gcd(pg,TV) = 1 and gcd(pr, 2) = 1. 

When < fc < TV, if 1 - e' 1 ^ = 0, then 2fc = TV, the 
real part of 



- fcrpir - fc(4g + r-)pn 

e 2JV -)- g JV 



is zero. 

If 1 - e* 15 *^ ^ 0, then 



■ krp-7, ■ fc(4g + r) P 7i 

e N + e N 



. k[4.q(N~l) + T]pT, 



. k{iq(N-l) + T]pw 

e « =0 



From equality (0, when < fci < TV, the inner product 
of the fcith row and itself of C(TV) is 1. We know that the 
product of C(N) and its transpose is an identity matrix. 
Thus, it is easy to get the inverse transform of transform ©. 
We omit it here. 



(2) 

We remark the following three points: 



(I) The transpose of C(N) is an orthogonal matrix, so the 
transform below is also orthogonal. 



inner product of the fcith row and fc^th row of C(N) is that 
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X(k) — a(n)x(n) cos 



71=0 



n(4:qk + r)pTT 
2N ' 
fc = 0,l,--- ,N-1. 



where gcd(pq, N) = 1, gcd(pr, 2) = 1, p, q and r are positive 
integers. 

(II) We can obtain a generalized DCT-IV-E transform, as 
follows: 



7V-1 



X(k) = a(k) x(n) cos 



n=0 



(2k + l)(4gn + r)p7r 
AN ' 
fc = 0,l,--- ,N-1. 



where gcd(pq, N) — 1, gcd(pr, 2) = 1, p, q and r are positive 
integers. 

(Ill) In a similar way, we can generalize discrete sine 
transforms, such as DST-II-E, DST-III-E, DST-IV-E. 
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— cos h cos 

2 27V - 1 27V - 1 

2(fci-fc 2 )37r 2(fei + fe 2 + 1)3tt 

4- cos • COS 
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Since e M = cos x + i sin x, we consider the real part of the 
following complex number: 



2JV — 1 -|_ p fc 2iV-l 



■ 2 fc( 2 JV-3)7i 

e 1 2 «-i 
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1 — e 2»-i 



III. A NEW TYPE OF DISCRETE COSINE TRANSFORM 

Let {x(n);n = 0, 1,2, N — 1} be a vector of real 
numbers. We define a new form of discrete cosine transform, 
as follows: 



X(k) 



27V 



-a(k) y * a(n)x(r- 



(2k + l)(2n + 1W 

■s , 

27V - 1 

fc = 0, 1, • • • . TV — 1, 



(4) 



When 0<fc<2JV-l,l - e 1 ^^ f 0. 
Note that 

, 4fc^(W-l) 
g fc 2JV-1 

2kn . , 2fc7r 

= cos(2K7r ) + z sin(2K7r ) 

V 2iV-r V 2iV + l ; 

. 2kir N . 2/C7T , 

= cos( 2iv3T)- iSm (^T) 

2A'ti 



By setting a = 2 n-i > we nave 



where 



o(jfe) 



1 else 

75 » = ^-l 



1 else 

The transform above can be written into matrix form: 

X(N) = C(N) ■ x(N) 
where X(N), x(N) are column vectors of length N, 



(2JV-3)ti 
2JV — 1 

1 

75 
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cos 2JV — 1 
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(2JV-3)(2jV-3j7i 
2JV-1 
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For the sake of simplicity, the coefficient 



27V-1 



" V2 

l 
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is omitted 



in the following discussions. Let < fci < k 2 < N — 1, the 



1 - e 2N ~ 1 

cos a + i sin a 

1 — cos 2a — i sin 2a 

cos a + i sin a 

2 sin a (sin a — i cos a) 

(cos a + i sin a) (sin a + i cos a) 
2 sin a (sin a — cos a) (sin a + zcosa) 
i 



2 sin a 
Thus, the real part of 
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From equality (0, for < k\ < &2 < N — 1, the inner 
product of the fcith row and foth row of C(N) is 



1,1 1. 1 

2(-2-2) + 2^° 



By setting a = ^N-i ' we nave 
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.. 2(2fc + l)(W-l)ir , 
(1 _ p 1 2AT-1 ) 
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2 2 sin a 

Now, for < ki < N — 1, the inner product of the fcith 
row and (AT — l)st row of C(N) is 

1, In / In, Is 

Thus, C(N) is an orthogonal matrix. Equivalently, trans- 
form dH) is orthogonal. 

It is easy to know that the product of C(N) and its 
transpose is an identity matrix. Thus, we can get the inverse 
transform of transform (jHJ. We omit it here. 

Similarly, we can obtain a new form of discrete sine 
transform as follows: 



X{k) = 



-N-l 



2N- 



n=0 



n) sin 



(2k + l){2n + 1)tt 
27V + 1 ' 
fc = 0,l,--- ,/V-l. 



And a new form of discrete sine-cosine transform as fol- 
lows: 

27V 



, v (2fc + l)(2n + l)7r 

( } = ^27vTT^ :c(n) cos " 
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fc = 0,l,--- ,AT-1; 
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n=0 



The orthogonality of above transforms follows from an 
analysis similar to that of new discrete cosine transform ([8]). 



IV. Generalized discrete W transform 

Bracewell (1983)advanced a discrete Hartlry transform 
(DHT) in the domain of real numbers. Its nuclear function 
cas wt = cos wt + sin wt is the sum of the real part and imag- 
inary part of Fourier transform nuclear function exp(iwt) = 
cos wt + i sin wt. Zhongde Wang (1984)generalized DHT and 
advanced discrete W transform (DWT), which is an unified 
form with two parameters. There are four meaningful types 
(including DHT) in DWT. The nuclear function of DWT is 
still cas wt = cos wt + sin wt. The discrete W transform has 



been used in frequency spectrum analysis, data compression, 
convolution computation and information hiding. Its theory 
and algorithms have received much attention for the last two 
decades (Wang,1989, Wang,1992, Kunz,2008, Zhang,1992). 

Let {x(n);n — 0,1,2, N — 1} be a vector of real 
numbers. The unified form of discrete W transform is defined 
(Wang, 1984, Wang, 1985) as follows: 

/yw-i 27r 



n=0 



AT- 

fc = 0,l,--- ,N-1 



The transform above is orthogonal only for (a, f3) £ 
{(0,0),(i,0),(0, It is discrete Hartley transform 

when (a,f3) — (0,0). The four transforms are defined as 
DWT-I, DWT-II, DWT-III and DWT-IV, respectively. 

Actually, the unified form of discrete W transform can be 
added with one more parameter 7, as follows: 

iyJV-1 

X ^ = VnT, < n ) sin ti + ( fc + a ^ n + 0)77]' w 



71=0 



fc = 0,l,--- ,N- 1 

1 r 



This is the proof that when (a, /3, 7) = (^,^,2pq) and 
gcd(pq, N) = 1, where p, q and r are positive integers, 
transform (|6]l is orthogonal. 

For the sake of simplicity, when (a, /3, 7) = (§, ^, 2pq) and 
gcd(pg, N) = 1, transform © is rewritten as follows: 



N-l 

X ^ = 7X7 E -^Wcas[(2fc + l)(qn + r)^], (7) 
v A n=Q 

fe = 0,l,--- 

where cas a; = cos a; + sin a;. 

Transform (O can be written into matrix form: 

X(N) = C(N)x{N) 
where X(N), x(N) are column vectors of length N, 
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where 77 = (2A^ — l)p7r. 

To prove that transform (O is orthogonal, or equivalently, 
to prove that H(N) is an orthogonal matrix. 

It is easy to know that 

casAcasB 
= (cos A + sin A) (cos B + sin_B) 

= cos A cos B + sin A sin B + sin A cos B + cos A sin B 
= cos(A - B) + sm(A + B) 



To simplify the problem, we omit the coefficient -j=. For 



I < ki < k2 < N, the inner product of the fcith row and fc2th 
row of H(N) is that 
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As for the function e lx — cos x + i sin a;, we consider 
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When < k < N and gcd(pq, N) = 1, we know that 
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hence 
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Thus 
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When < k < 2N, if 1 - e 4 ^^ = 0, then k = N, thus, 
the imaginary part of 



. 2k r pw . 2fc(g->>7t 
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,- 2fc[ g (JV-l) + r]p7r 



is zero. 



When < k < 2N and if 1 - e l ^~ ± 0, then 
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Thus, 
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From equality (0, the inner product of the fcith row and 
feth row of H(N) is zero for 1 < hi < < N, namely 
H(N) is an orthogonal matrix. 



For 1 < fci < N, the inner product of the fcith row of 
H(N) and itself is that 

^{[cas( 2fcl ZpIElf + [cm (2fci-l)( g + rW ia + 



AT 
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Note that < 2k 1 - 1 < 2JV, from the proof above, we 
know that 

. 2(2fci-l)r P 7r 2(2fci - l)(q + r)m 
sm — h sm 



N N 
2(2fci-l) [g(JV - 1) + r]p7r 
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= 



It follows that the inner product of the fc x th row of H(N) 
and itself is 1 for 1 < hi < N, namely the product of H(N) 
and its transpose is an identity matrix. Thus, it is easy to get 
the inverse transform of transform (0, 



x(k) 



N-l 

-±= £x(n)cas[(2n + l)( g fc + r)^], 

* n— 



fc = 0,l,--- ,N-1, 



where cas x — cosx + sinz. 

We should remark the following two points: 

(I) There is no limit to parameter r in the new transform 
(2), r is an arbitrary positive integer. 

(II) In a similar way, we can generalize DWT-IV as follows: 

1 JV_1 

X ^ = 7^ E *Wcas[(2fc + l)(2qn + r)& 

fc = 0,l,--- ,JV-1 

where gcd(pt7, iV) = 1, p, q and r are positive integers. It is 
degraded to transform (2) when r is an even number. 

Consider two special cases of transform ©: 
When (a, /3, 7) = (5,5,4) and iV is an odd number, 
transform (|6]l can be written as follows: 

y ' ./at V ' AT ' 



n=0 



TV 

fc = 0,1,--- ,N- 1 



When (a, /3,7) = (5, i , 2g) and gcd(q,N) = 1, transform 
© can be written as follows: 

jV-1 



(2fc+ l)(qn+ 1)tt 
2iV ' 
fc = 0,l,--- ,iV-l 
We know that transforms above are orthogonal. 



X(k) = —= a;(n)cas 

v n— 
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